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Abstract 

We present a general class of noncolinear colliding wave solutions of the 
Einstein-Maxwell equations given in terms of fourth order polynomials, which 
in turn can be expressed through Jacobi functions depending on generalized 
advanced and retarded time coordinates. The solutions are characterized by 
six free parameters. The parameters can be chosen in such a way to avoid the 
generic focusing singularity. 
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I. INTRODUCTION 



Recently, the collision of plane-fronted gravitational waves possibly coupled with elec- 
tromagnetic waves, has been extensively studied Because gravity is always attractive, 
it was expected that focusing of the waves would occur and one of the interesting questions 
is how much focusing does general relativity predict. Within this framework, strong focusing 
would appear by the development of spacetime curvature singularities. Many solutions has 
been presented so far, describing collisions of plane-fronted gravitational and electromag- 
netic waves. And quite a few of them do develop Cauchy horizons. It is important to stress 
the fact that all of these solutions contain only second degree polynomials and that till now 
nobody has gone ahead to higher degree polynomials in the sense we are going to explain 
below. 

All known colliding wave solutions are characterized, in null coordinates, by the line 
element 



ds = e dudv + gabdx"'dx , a, 6 = 1,2, (1-1) 

where x^, are ignorable coordinates and metric functions depending on {u,v) 0. When 
written in [rj, fi) coordinates 0, where rj is a. measure of the time from the instant of collision 
and /i is a measure of the distance between the plane wavefronts, the metric becomes 

ds' = e'^iv' - /i') - ^) + 9a,dx'^dx\ (1.2) 

where the polynomials P{ri), Q{n) are at most of the second degree in their respective vari- 
ables. Generating techniques and other methods have allowed a great deal of generalizations 
in relation with the metric components gab (sometimes called Killingian sector) HBHTSI. The 



generated solutions preserve the seed's structure of the polynomials P{rj), Q{fi), changing in 
general the form of the metric function 7. 

In this paper, we pursue further and present solutions with P{ri) and Q(/i) being polyno- 
mials of fourth degree, by means of a generalization of the advanced and retarded coordinates 
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concept. This generalization endows the corresponding spacetimes with a geometric struc- 
ture for the time and independently for the spatial coordinate richer than the already known 
colliding wave solutions. In the limit of polynomials of second degree one recovers the widely 
used retarded and advanced time coordinates {u,v). The main purporse of this paper is to 
give the formulation of the most general class nowadays known of colliding waves with fourth 
degree polynomials. As usual, it is assumed that in the corresponding spacetime, the two 
waves approach each other from opposite sides in flat Minkowski background; after the col- 
lision, a new gravitational field evolves, which satisfies certain continuity conditions. The 
colliding plane waves possess five symmetries, while the geometry resulting after the collision 
has two spacelike Killing vectors. 

The plan of the paper is as follows: In Sec. II we revisited some generalities on colliding 
waves concept. In Sec. Ill the general representation through advanced and retarded time 
coordinates for fourth degree polynomials is introduced. In Sec. IV the class of colliding 
waves with fourth degree polynomials is presented. In Sec. V the results are discussed. 



As was stated previously, in this section we review briefly some generalities on electrovac- 
uum colliding waves, introduced by Ernst et al. p|. 

The set of colliding waves solutions is described by the line element 



where = x, x"^ = y are ignorable coordinates. The domain of the coordinate charts consists 
of {x,y) G and {u,v) G R^; it is the union of four continuous regions: I := {{u,v) : < 
M < 1,0 < < 1} , // := {{u,v) : u < 0,0 < V < 1}, III := {{u,v) : < m < 1, t; < 0}, 
IV := {{u,v) ■.u<0,v < 0}. 

In the region IV, a closed subregion of the Minkowski space, it is required that 



II. COLLIDING WAVES 




a,b= 1,2, 
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where is the electromagnetic vector potentiaL By scahng-shifting transformations 
the metric can be brought to standard Minkowski metric. In region //, the met- 
ric components and the electromagnetic vector potential depend only on v i.e. g^^p = 
gfj.u{0,v), = y4^(0,f). In region III these fields are functions of the coordinate u, i.e. 
9fiu = QfMuiu, 0),A^ = Af^{u, 0). In region I, which is occupied by the scattered null fields, the 
metric components and the electromagnetic field are functions of both u and v coordinates. 
In this way, the problem is reduced to know the metric (coframe) and the electromagnetic 
vector fuction A^. 

A. Newman-Penrose tetrad 

We introduce a complex null frame according to the conventional description: 

e„ = (1, n, m, m) (2.3) 

with 

1^ = = = = lm = nm = lm = nm = 0, (2-4) 



1 ■ n = 1, m-m = — 1. (2.5) 

The basis of the corresponding 1-forms, the coframe, is denoted by i?". If we lower the 
coframe index, then we use the same letters as with the frame: -^a = (1, n, m, m). Then the 
coframe reads 

= (n, 1, -m, -m) . (2.6) 
According to ( p.4[ ) and (p3|), we have 

g = 2(^^^ (^^^ -^^ (g)^^) =n^p^'' ^^f^ , (2.7) 

with 
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'^O 1 ^ 



10 
-1 
0-10 



(2.8) 



B. Self dual part of the complex Weyl 2 form 

Quite independently from the type of coframe we are using, we can consider some prop- 
erties of the Weyl 2-form. Let us denote the Weyl curvature 2-form by 

C^p = ^C,,^p^^Ar. (2.9) 

If, in four dimensions, an arbitrary complex 2-form u is given, then its Hodge dual *u 
is again a 2-form. Therefore we can build linear combinations of both: 

LU^ ■=^{u±i*u) . (2.10) 

Here i is the imaginary unit. We call a;"*" the self-dual and uj~ the antiself-dual part of the 
2-form uj. Clearly, 

uj = uj+ + u}~. (2.11) 

A form UJ is called self-dual, if its antiself-dual piece vanishes, i.e., if u; = it is called 
antiself-dual, if a; = —i*uj. 

Let us now assume that the form uj is real, that is, = cJ, where the bar denotes complex 
conjugation. Then, according to ( |2.10| ), they fulfill the relations: 



uj^ = uo , and uo~ = uj^ . (2.12) 

Consequently, a real 2-form uj = uj can be alternatively encoded into its self-dual part uj^. 
Incidentally, the same is true for its antiself-dual part. 

The Weyl 2-form is a real quantity. Accordingly, we can also take its self-dual part 



C^p--=\{C^p + i''C^p) . (2.13) 



Expanding the Weyl 2-form C^^ in terms of the basis ■t?'^ A we have 



ctp = \cUp^'^^\ (2.14) 



which we will asumme to be in the standard 2-form selfdual basis [U, V, W]. The most gen- 
eral form of our Weyl field components, compatible with colliding wave spacetime structure 
is given by |p!4| : 



= A^{u, v), Ct^^ = 2^oUabUcd + 2^2{UabVcd + VabUcd + WabW^d) + 2^ ^VabVcd , rcgiou I, 
= A^{v) , C+^a = 2^oUabU,d , region II , 
A^ = A^{u) , = 2^,VatV,d , region III , (2.15) 



with 



Wab = maifib - mbfha - kJb + hL , 

Vab = kaTTlb - huia , 

Uab = -Lrhb + hrfia ■ (2.16) 
where rria, m^, ka and la are null tetrads. 

III. GENERALIZED ADVANCED AND RETARDED TIME COORDINATES 

We shall consider a metric of the form 

'-^-'iwrm)^ <'-^' 

where P, Q are fourth degree polynomials in their respective variable, and 7 depending on 
both (p, q) coordinates. 

Since we have real polynomials of fourth degree, when looking for advanced and retarded 
time variables one has to deal with elliptic integrals, i.e. we shall use the Legendre first kind 
integrals: 



dv 

/ : . = /^^(0' k) , (3.2) 



where the fourth degree polynomial 6*4 (r), is represented by the product of four monomials 
containing the corresponding roots, the coefficient of the higher degree (fourth), denoted by 
04, is assumed to be 04 = ±1; if this were no the case, i. e. 04 = ±04, 0:4 > 0, then dividing 
by 04 one arrives at the above equation where now fi has to be interpreted as in what 
follows we shall adopt these conventions when needed. The real roots are denoted by rj, 
with j = 1, . . . , 4 and ri > r2 > > r^, while the complex roots by Si ± iti, S2 ± it2, with 
Si > S2, ti > t2 > 0, we introduce also the following useful notation 



(r,/3,7,5) 



rik = rk-ri, (z, /c = 1, 2, 3, 4) 
r — 7 P — S 



r — 6 /9 — 7 ' 



tan Oi = , tan Uo = , 

^1 ~l~ ^2 /I ^1 — ^2 
tan ^3 = , tan 64 = , 



tan 



Si - S2 Si - S2 

COS 6*3 

COS 6^4 ' 

z/ = tan [(^2 - 0i)/2] tan [(^2 + ^i)/2] . (3.3) 



The elliptic integral F{(j), k) is the Legendre integral of the first kind. As it is well known, 
the standard form for writing this function is 

■■F{(j),k), (j) = 8imz, (3.4) 



1 — k"^ sin^ (, 

where amz denotes the function amplitude of z. Replacing through u according to 

cj = sinf/) = sin(amz) := sn^; , z= , = = Fju, k) , (3.5) 

^0 ^(l-cu2)(l-Fa;2) 

where snz belongs to the Jacobi family of elliptic functions (sn^;, cn 2;, dn z), with well es- 
tablished analytical properties. The values and main properties of these function can be 
found in [p!5| j. 

Going back to the line element ( |3.1| ) and assuming that P{p) and Q{q) are polynomials 
up to fourth degree on p and q respectively, then the two dimensional line element, first two 
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terms in can be written in terms of the retarded and advanced time coordinates u and 

V, respectively, as follows 



-27(P.9) 



\W)~qW)I 



-^^""^"^dudv. 



(3.6) 



with 



du 

dv 



dp ^ dq 



± 



de 



P4{p) \fQM v^^/r^Esin^ V^Vl-fclsin^^ 



(3.7) 



where fj.(j,, are constants and 04, /?4 are the moduli of the coefficients of the higher power 
in the polynomials P and Q, correspondingly. The variables u and v range the values from 
—00 to +00. 

The integration of ([3171) yields 



u 



(3.8) 



where — cxd<m<oo,— oo<^;<c)0. Hence, one has relations between the Legendre elliptic 
integrals and the null coordinates u, v, i.e. 



F{<P,k^ 



u + v) 



F{e,kg 



[U — V] 



(3.9) 



Consequently, the inversion formulas bring the following functional dependence upon the 
generalized advanced and retarded time coordinates 



am 



9 = am 



2/i0 



u + v) 



[u — V 



sn 



Q = sin 6 = sn 



2/i,/, 
2^e 



u + V) 



[u — V] 



(3.10) 
(3.11) 



which allow us to give the original p and q coordinates, initially expressed through and 9, 
in terms of u and v null coordinates, i.e. 



p = p{u,v), and q = q{u,v) 



(3.12) 



Moreover, one can introduce a new set of variables {u, v) 



u = sin y—j , and v = sin ^-j . (3.13) 

where the ranges of the variables are — oo < -u < oo, — oo < -y < oo and — 1 < m < 1, — 1 < 
f < 1. These variables are such that 

du dv , ^ 

du = 2—, dv = 2y (3.14) 

and 



arcsinw ± arcsinf = arcsin[Ml/ ± vU] (3.15) 



2 

where U = y/l — m^, V = y/l — and the corresponding term in the line element 

7~ 7~ Adu dv ,^ 
du^dv = 4— ^ y. (3.16) 

Having in mind colliding wave interpretations of the results, one restricts u and v to the range 

0<'U<l,0<f<ltobe closed to the standard colliding wave approach. Depending on the 

character, real or complex, of the roots of the fourth degree polynomials, one may encounter 

in general combinations of different possible cases. It is always possible to introduce the 

retarded and advanced time coordinates u and v. However, only certain solutions satisfy 

the Ernst ^ requirements for colliding waves. We shall restrict ourselves to metrics which 

can support a colliding wave interpretation. 

It may occur several mixed cases depending on the degree of the polynomials appearing 

in the non-killingian sector, for instance, 

d]^ _ d(f dj^ _ dq^ ^^^^^ 
P4 Q2 -P3 Q2 

where the subindice indicates the degree of the polynomial we are dealing with; the possibility 
of interchanging P by Q and p by g and viceversa has to be taken also into account. The 
case ^ ~ ^ can be treated in the manner we did with the fourth degree polynomial. 



Let us consider the case of — Introducing the null variables 



^2 



±^^, (3.18) 
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in the case when Q2 can be brought to the form 1 — q^, they become 

► = ii4>F{(j), k^) ± arcsing, , 



u 

V 



therefore 



u — V 

_ 2 ^ 
u-\-v 



arcsmg, 



/'t' dx 



sm^ X 



fi^z, (j) — amz, 



hence 



u — V U + V 

g = sm(^— ), p = p(^— ) 

where the relation of p in terms of u, v depend on the roots of the polynomial P4. 
case with P3, Q2 can be treated in similar way. 

For completeness we also briefly comment the following cases 



dp'^ dq^ dp^ 



PaAp) QM p^M 

the corresponding v are 



-(3{dq-^)\ g>0,/3>0, 



u 

V 



> = ii^F{(p,k^)±i3q^, 



thus 



(ii) 



,u + v. .u-v.o 



dp"^ dq^ dp^ dq^ 



P2{p) Q2{q) 



the corresponding u, v are 
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u 

V 



arcsinp ± arcsin q, (3.27) 



hence 



p = sm(^), g = sin(^), (3.28) 



(iii) 



dp 



2 



the corresponding v are 



dq\ (3.29) 



u 

V 



> = fJ.^F{(j), k^) ± g, (3.30) 



hence 



2 ^' ^ 2 ' 

(iv) 



P = P{^^): (1 = ^:;—: (3.31) 



-/5(d,^)^ /5>0,g>0, (3.32) 



P2{p) Qi 

the corresponding u, u arc 



arcsinp ± ^q^ , (3.33) 



hence 



(v) 



Pi{p) QM 

the corresponding S, v are 



p = sin(^— ), g=(^^), (3.34) 



"^^^ "^^^ (dp^)^-(^ig^)^ p>0,g>0, (3.35) 
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(3.36) 



thus 



. n ) ' 



,U 



(3.37) 



A. Different transformations of p[u,v) and q{u,v). 

In this subsection, we present all possible cases depending upon the roots of the polyno- 
mials Pip) and Q{q) and the corresponding transformation. 

1. For P and Q with all real and different roots, i.e. 

P = Q;4e± [p - pi] [p - P2] [p - Ps] b - P4] , e± = ±l, q;4>0, (3.38) 
Q = 13av± [q - qi] [q - ^2] [q - ^3] [q - ^4] , i/± = ±1, /34 > , (3.39) 

The corresponding transformations read as follows 





^± 


Inter\'al 


= 


A-2 


four 
real 
roots 


e+ 


Pi < p or p < p4 


P1P42-P2P41 ^Sn2^(«+-u) 


2 


(«!, 0:2, a4, as) 


P42-P4lfsn^(u+-!;) 


2 


e+ 


Pa < P < P2 


P3P42 -P4P32 (™ 2^ 


2 


P42 -P32 ( '^'^ (w+f ) 


2 


e_ 


P4 < P < P3 


P4P31+P1P43 (^sn2^(u+l>) 


2 


(as, q;2, "4, tti) 


P31+P43 fsn^ («+■!;) 


2 


e_ 


P2 < P < Pi 


P2P31 -P3P21 (^"^ 2^ {u+v) 


2 


P31-P21 ^sn 5^ j 



In the transformations given above we must understand that u — 2 arcsin u and v — 
2 arcsin V, such that u±v — 2 arcsin(rty±^;C/) and //^ has to be understood as //^/y^ 
when 0:4 7^ 1. For Q{q) the same Table is valid making the changes: p ^ q, Pi q^, 
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e± — > i/±, ijl^ — > ijL0 and + ^; — > w — t), Si ^ ai and — > Tj. The same applies in the 
following cases. 

2. For the case in which P and Q possess two real and different roots, and two complex 
roots, i.e. 

P = aie±\p-pi]\p-p2]\p-{si + iti)]\p-{si-iti)], 04 > 0, (3.40) 

Q = fe[g-gi][?-?2][g-(?i + iTi)][?-(ai-ifi)], /34 > 0, (3.41) 

the transformations are given by 



Pa{p) 



Interval 



p{u,v) 



two complex and 
two real roots 



-00 < p < 00 



P1+P2 
2 



ti/e-cnirj—(u+v) 



(sin(^))' 



3. In the case P and Q with four complex roots, respectively, i.e. 

P = q;4 [p - (si + iti)] [p - (si - i^i)] b - (S2 + ?i2)] [p - (S2 - ii2)] , (3.42) 

Q = (3A[q- {ai + iri)] [q - (ai - iri)] [q - {02 + ir2)] [q - (^2 - ^72)] , (3.43) 

the transformations are now 





e± 


Interval 


p(w, v) = 




Si > S2 


e+ 


—00 < p < 00 


Si + ti 


'^'^2^ («+t')cos ( ^ + ^ ) +cn («+«)sin ( ^ + ^ ) 


sin^ 6*5 


_cn^ (^+^;)cos( ^ + ^) -sn^ («+^)sin( ^ + ^) _ 


Si = S2 
Bti>t2 


Si - ti 


^-2^(^+^)-°Kf + ^)+--2|^(^+")«-(^ + ^) 




_^°2^(^+")^°<^ + ^)-«°2^(^+">i°(^ + %)_ 



We include the transformations for the case of polynomials of third degree: 
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e± 


Interval 


p{u,v) = 




three 


e+ 


P3<P<P2 


P32 
P3 




P32 
P31 


real 




Pi<P 


Pi 








1 






roots 


e_ 


P<Pi 


Pi 1 




P21 




P31 






P2<P<Pl 


P2P:n- 


-P3P21 (^sn^(«+D)^ 








P31- 






two 


e+ 


Pi<P 


Pl 






COS $1 ^1— cn^^ {u+v)^ 




complex and 










(sin 


'01 I tt' 




one real 


e_ 


P<Pi 









For the polynomial Q the same Tables are valid making the changes: p ^ q, Pi ^ 
qi,e± — > z/±, /i0 — i> fig and S + ^ -u — -i;, Sj ^ cij and —>■ Ti. 

All possible combinations between the transformations corresponding to the different 
possible ranges of p and q should be taken into account. Moreover, it is straightforward to 
perform such combinations by using the explicit transformations in tables given above. 

B. Limiting transition 

For completeness we include the limiting procedure to obtain from the fourth order 
polynomials the third and second degree ones. 

P = e±a^{p - pi){p - p2){p - p^){p - Pa) 
P 

= ami l)(p-P2)(p-P3)(p-P4) 

Pi 

p 

= a^i l){p - P2){p - P3){P - Pa), (3.44) 

Pi 

Taking the limit when pi — > cxd we recover a third degree polynomial: 

lim P4 = e±a4{p - p2){p - p^){p - Pi) = e^a^Ps, (3.45) 
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Analogously we take the limit on P3 when p2 — ^ 00 to get a second order polynomial: 
lim P3 = lim I £±04^2 (— - - P3){p - Pa) \ 

= lim \e±a3{— -I){p-p3){p-p4)\ = e±a3{p-p3){p-p4) = e±a3P2, (3.46) 

P2-»0O \ P2 J 

The corresponding procedure can be achieved for Q4 to obtain a second order polynomial. 



C. Concrete example 

As a concrete example let us consider, for instance, the case P with four real and different 
roots and Q also with four real and different roots, in which the coefficients of the higher 
degree are 1, and the roots of P and Q are denoted by Pi and respectively. For the elliptic 
integral depending on the p-coordinate, with P(p) = {p—Pi){p—P2){p—P3){p—P4), (e± — !)> 
one has the relation 



/P{p) 

where the explicit relation between p and 4> reads 

P1P42 - P2P41 sin^ (f) 

• 2, ^ 

P42 - P41 sm 

. 2 , P42 (P-Pl\ , . ( , PA2 (P-Pl 

sm = — , — > = arcsm ±. — 

P41 \P-P2 ) \ \PAI \P-P2, 

The parameter A;,^, < A;| < 1, is given by 



(3.48) 



,2 / N Pi-Pa\ P2-P3\ 1 , . 

^ = Pi,P2,P4,P3 = , and = — ^ = ^ . 3.49) 

\Pl-P3J \P2-P4J VP31P42 

It is straightforward to find the analogous expression for the elliptic integral depending on 
the g-coordinate. The explicit expression of the coordinates p and q through u and v is the 
following 

~^ P1P42-P2P41 {Sn^ju + V))' 

p{u, v) = — —2— , (3.50) 

P42 -P41 [snj^{u + v)) 

542 - 541 (sn2^(M-'u)j 
16 



In terms of the null variables u and v which bring the line element to the form 

dp^ dq^ dudv 

-TT- = dudv = 4— — — 

P4 Q4 uv 



the p and q variables amount to 

P1P42 — P2PAI ( sn— (arcsinfuV + vU)) 



p{u,v) 



q{u,v) 



Pa2 — Pil (sn2^(arcsin('uV^ + vU)) 



2 ' 



gig42 - g2g4i (sn2^(arcsin(uF - vU)) 



2 ' 



^42 — qAi (^sn2^(arcsin(My — vU) 
On the other hand, the relations between and 6 with p and q are 



±6 = arcsin , 



P42 P- Pi 



± arcsin 



g42 q- qi 



(3.52) 



(3.53) 
(3.54) 



(3.55) 



\PA1\P-P2j '\qii\q-q2j' 
where and 9 stand correspondingly for the arguments of the elliptical Legendre integrals 
related with the integration of P and Q. We define the auxiliary non-null variables u' and 



V as 



u = sm 



-9^ 



V = sm 



(3.56) 



with 



arcsin m' ± arcsin t;' = arcsin u'vl — v^ ± v'vL — u 
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arcsin \iiV' ± v'\J'^ 



(3.57) 



where 



^' = Vl - v''^ , and f/' = Vl - 



(3.58) 



Thus, when we substitute ( p. 571 ) and ( |3.58| ) into ( |3.5(]| ) and ( p.51| ) it is staightforward to find 



p(m', V 



q{u,v') 



P1P42 - P2P41 {u'V + v'U'Y 
P42 - Pai (u'V + v'U' 

qiqi2 - 92^41 {u'V - v'U') 
q42 - q4i {u'V - v'U'f 



(3.59) 
(3.60) 
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Moreover, from ( pl7|) and ( p.57|) , the advanced and retarded time coordinates u and v are 
related with the variables u' and v' through 



du 
dv 



± 



+ 



yJl-kliu'V + vU'y ^l-kl{u'V' ~v'U') 
1^4) _^ l^e 



du' 
IF 

dv' 
V' 



(3.61) 



'1 - kliu'V + v'U' f ^1 - kl{u'V' - v'U')\ 

According to our general procedure ( p.61| ) can be inverted, yielding u' = u'{u,v) and 
v' = v'{u, v). 

In the limit of polynomials of second degree , i.e. P2, Q2, one has k^ = kg = 0, and 
fJ'ij), fJ^e — ^ 1, hence one recovers from ( p.61|) the widely used retarded and advanced time 
coordinates 

du' dv' 
du = 2— dv = 2—, (3.62) 

In the next Section we present an explicit class of colliding wave solutions. 

IV. A CLASS OF COLLIDING WAVES 

Let us consider an explicit solution of the Einstein-Maxwell equations in region I, given 
by the following coframe in coordinates {x,p, q, y): 

^' = j^^dq, ^^ = j^^idx + p'dy), 



(4.1) 



Here we have the functions H = H{p,q), P = P{p), Q = Q{q), and A = A{p,q). The 
coframe is assumed to be orthonormal 



(4.2) 



Then the metric explicitly reads 
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l_i^^(^dx + p^dyy dq^ + ^dp^ + ^ (dx - q^dyf^ . (4.3) 



with 



H{p,q) := 1 -/ipg 

P{j)) := h — + 2np — ep^ + 2m^p^ + \ ~^ ~ l^^b — jJ^e^] p^ 



Q{q) ■= -{h + e^) + 2mq - eq^ + 2n/i + (^^ + /"^^ " /"^^^ j ^ > 
Ab,g) :=p' + g', (4.4) 

while the nonvanishing electromagnetic field components are given by 

_ {p^ -q^)e- 2gpq 

F ■= v^F 

_ g{q^ -p^) -2epq 

J^xp ■ ^2 

Fyp := -q^F^p . (4.5) 

In order to express our solutions ( |4.4| ) and ( |4.5| ) in terms of the advanced and retarded time 
coordinates u, v, one has to use the inversion procedure,presented in the previous section, 
yielding p = p{u,v) and q = q{u,v). The explicit representation requires the use of tables 
for the sn function. 

Our class of solutions, defined in region I, can be extended to the full spacetime by 
introducing the Heaviside step function 

r 1, X > 

Q{x) = , (4.6) 

[ 0, X < 

with 02 (x) = 0(x), and performing the following transformation 

u—^Q{u)u, and v Q{v)v . (4.7) 

For the case of real roots of the polynomials, treated extensively in Sec. III.B, we obtain for 
region II 
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p{u,v) 
q{u,v) 



Pm2 - P2P41 



P42 - P41 {s^2^v) 

g42 -g4i (sn2^?))^ 



(4.8) 
(4.9) 



where 



sn— = sn 
2/i 



arcsm v 



(4.10) 



For region III 



p{u,v) 
q{u,v) 



Pm2 - P2P41 (^sn^^M 

gig42 -g2g4i (sn^^w)' 

g42 -g4i (siis^u)^ 



(4.11) 
(4.12) 



where 



u 

sn — = sn 
2/x 



■ arcsm m 



(4.13) 



At this point it is important to mention that not every cyhndrically symmetric spacetime, 
even in vacuum, satifies the requirements of Ernst |^ for being interpreted as colhding 
waves, only certain classes of cyhndrically symmetric solutions can be thought of as solutions 
generated by collision of waves. 

Moreover, up to here our analysis has been concerned mostly with metrical aspects 
leaving aside the task of having suitable field and matter energy-momentum tensors. As 
it is well known, each participating wave in the head-on collision is assumed to be plane 
fronted, i.e. it is characterized by a covariantly constant null eigenvector k, i.e., kf^.^, = 0, of 
the Weyl Petrov type N conformal tensor. This condition implies a null energy-momentum 
tensor for the electromagnetic field, i.e. a radiation field or electrovacuum, — let us call it 



null field for short — into the Einstein field equations, c.f. Kramer et. al. [0], Section §21.5. 
Therefore, the existence of colliding wave solutions with non-vanishing cosmological A term 
is thus forbidden. 
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Usually one derives colliding wave solutions using the backward way procedure, namely, 
one starts with a cylindrically symmetric solution and an energy-momentum tensor in terms 
of fields depending on the usual two null variables, say u and v, restricted to a closed interac- 
tion region, region I in our approach. Furthermore, by accomplishing on it the transforma- 
tions ( [4. 7] ) with Heaviside functions one obtains the spacetimes prior to the collision. These 
spacetimes occur to be type N solutions to the Einstein equations in vacuum or coupled with 
null fields, regions II and III in our approach, where the gravitational and null fields depend 
solely on a single variable {u or v). At this level, the colliding wave condition allows only 
for fields which can be matched to the fiat Minkowski spacetime (region IV) background 
where the wavefronts propagate [|. In other words, if the cylindrically symmetric spacetime 
is endowed with a cosmological A term, one has to set it equal to zero in order to connect 
smoothly the different regions of the collision and in order to avoid meaningless jumps of 
the cosmological constant across the null lines. 

On the other hand, it is worthwhile to mention the results by Chandrasekhar and Xan- 
thopoulos |]T^ concerning colliding wave solutions with a perfect fluid (stiff matter, fluid 



pressure = fluid energy density). They found that after the collision takes place, the incom- 
ing plane waves supporting radiative flelds, give rise to solutions with stiff matter perfect 
fluid tensor. It is an example of the transformation of massless particles describing null 
trajectories into a perfect fluid in which the stream lines describe timelike trajectories. 
In full, in order to interpret our class of cylindrically symmetric solutions ( [4.4| ) as a 



colliding waves class of solutions, it is compulsory to set A = on it [16 



V. DISCUSSION 

In this paper we present the formulation of the most general class nowadays of colliding 
waves with fourth degree polynomials determining the non-killingian sector, by means of a 



^We appreciate the referee's comment which draw our attention to this respect. 
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generalization of the advanced and retarded coordinates concept. As it has been pointed 
out, after setting the cosmological A term equal to zero, this class of solutions describes the 
scattering of two noncolinear polarized gravitation plane waves, having a null electromag- 
netic field as a source. At the leading edge of each colliding type-N gravitational wave, 
the curvature tensor exhibits delta and jump discontinuities. The former is interpreted as a 
gravitational impulsive wave, whereas the latter is attributed to a gravitational shock wave. 

As stated above, this class of solutions, defined again in region I, can also be extended 
to the full spacetime by introducing the Heaviside step function 6. 

The electromagnetic field present delta singularities and jump discontinuities. However, 



the Bianchi identities hold in a distributional sense, see [18 



There are no problems on the right-hand side because the delta type singularities of 
the curvature are multiplied by the smooth distributions \Jl — Q{u)u^ and y^l — 6(f )f ^, 
respectively. 

On the other hand, with respect to the singularities of the solution, we have in Region I, 
the solution is of type D in the classification of Petrov, this is so because the only novanishing 
Weyl scalar is \l/2, and explicitly, it is given by 

+q{q^ - 3p^)(a3 - mi) - - 3g^)(ai + mg)] , (5.1) 

where the parameters ai = 2n, = 2m. The parameters e and b do not appear in the 
curvature scalar, however the properties of the solution depend on them, as we shall see. 
The complex structure of \l/2 is due to the noncolinear polarization of the incoming waves 
(nondiagonal metric, see |]19[)- Since we have the explicit expression for the Weyl scalar, 
Eq. ( |5.1| ), we can study the behavior of the invariants of the spacetime and determine if the 
solution develops a curvature singularity a finite time after the instant of collision. Penrose 
attributed the development of the singularity in most colliding wave solutions to the self- 
focusing of the gravitational waves on scattering. In this case the invariants are given by 
I = 3<ifl and J = ^l- Then the singularities of \l/2 are clearly the same for the invariants. A 
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singularity occurs when A = + = 0, i. e. when simultaneously p = and q = 0. The 
interesting question is if there is any case in which the nature of the roots do not permit 
that simultaneously p = 0,q = 0. If this is the case, then the solution is free of the focusing 
singularity. A brief analysis shows that this occurrence can be avoided for certain values of 
the parameters of the solution. 

Let us consider the expressions given in Sec III.C for the case of P and Q both with four 
real roots. In this case (p.53|) and (|3.54|) express p and q in terms of {u, v). Analyzing those 
expressions we see that p = 0, q = simultaneously if, for instance, p42 = 0,p4i = and 
Q42 = 0, g4i = 0, i. e. when at least three of the four roots coincide: p^ = P2 = Pi and also 
q4 = q2 = qi- The assumption of the occurrence of three equal roots for P and for Q leads to 
algebraic conditions that must be satisfied by the parameters of the solution, b, g,n, 171,6,6 
and fi. Some of these conditions are incompatible. For instance, in the case considered, the 
following condition must be valid: + e^) = — g"^). This is only true in the case if 
+ (7^ = 0; i. e. if + (7^ 7^ we do not have that p = and q = simultaneously and 
then the singularity fails to occur. Then we conclude that in some cases the existence of the 
electromagnetic field permits to tune the parameters in such a way to avoid the focusing 
singularity. A related result was obtained already by Chandrasekhar for a colliding wave 



solution with a scalar field ||2^ and by Garcia et al. in the framework of metric-affine gauge 



theories of gravity [21,22 
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